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We present a leading-order pQCD calculation of the helicity-fhp A — > N-f* matrix element 
Go (Coulomb quadrupole amplitude C2), taking into account the transverse momenta of the 
quarks and the contribution from the gluons. In the large Q 2 limit, its scaling behavior 
acquires a double-logarithmic correction log 2 (Q 2 /A 2 ) compared with the standard scaling 
analysis, due to the contribution from the orbital motion of the small- a; partons. Based on 
this and on the latest JLab experimental results of the C2 — Ml ratio Rsm at Q 2 — 3 ~ 4 
GeV 2 , we make a phenomenological prediction for the latter at higher values of Q 2 . 



In recent years, there have been continuing interests in the electromagnetic N — > A transition. 
One of the earlier interests was due to the work of Becchi and Morpurgo |lj . In that work they had 
shown that in the context of the symmetric, non-relativistic, SU(6) quark model, the transition 
N — » A is a pure magnetic dipole Ml and the contribution from the electric quadrupole E2 is 
zero. One crucial assumption in their derivation of this 'selection rule' is that the quarks in both 
the nucleon and the delta are in the zero orbital angular momentum states. 

These predictions were to be considered as a check to the validity of the quark model, which 
was still questionable in those days. Later experimental measurements showed that, indeed, the 
magnetic dipole Ml contributes predominantly to the N — > A transition, while the contribution 
from E2 is small but non- vanishing. 

The non-vanishing value of E2, and also the Coulomb quadrupole C2 in the case of a virtual 
photon, has generated much theoretical interest. One way to account for E2 ^ is through 
the D-wave mixtures in the ./V and A wave functions 3]. Another way is through the two-body 
electromagnetic currents from one-gluon and/or one-pion exchange between constituent quarks [J]. 
In the latter case, it is argued that E2 transition is, mainly, due to a two-quark spin-flip operator. 
On the other hand, in the large N c limit of quantum chromodynamics it has been shown ||| that 
Rem = E2/M1 is of order 1/N 2 . To derive this result, no assumption about orbital angular 
momentum of the quarks was necessary. More recent work in this direction can be found in Ref. 



Another major issue related to the N — * A transition, and in general to any hadronic exclusive 
process, is the applicability of perturbative QCD (pQCD) at the range of values of momentum 
transfer Q 2 accessible in the current generation of experiments. In terms of the ratios Rem and 
Rsm = C2/M1, pQCD power counting predicts that 0], in the limit Q 2 —> oo, Rem —> 1 and 
Rsm — ► const., up to logarithmic corrections to be discussed in this paper. The former prediction 
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has not yet been observed experimentally. In fact, up to Q 2 = 4 GeV 2 , Rem stays negative and 
very close to zero 0,0,0,0,0]. The comparison between data and pQCD prediction for Rsm 
is the main topic of this paper. 

To make the pQCD predictions more relevant at finite Q 2 where data have been and will be 
taken, one has to go beyond the asymptotic power counting, make detailed pQCD calculations of 
hard scattering amplitudes and derive the factorization formula for the experimental observables. 
In the present case, the relevant quantities are the three independent helicity matrix elements [3]: 

Ga = ^-(A,A a |£(A)- J\P,X p = 1/2) , (1) 

where Mn is the mass of the nucleoli, A = 0, ±1, and Aa + 1/2 = A. A pQCD calculation of 
the helicity-conserving matrix element G + has been known for many years However, a 

pQCD calculation of the helicity non-conserving amplitudes has not been explored in the literature 
until recently, because it necessarily involves the orbital motion of the quarks 

0,0- 

In this paper, we perform a pQCD calculation for the helicity-flip matrix element Gq. The 
technical details are basicall y si milar to those in the calculation of the Pauli form factor ^(Q 2 ) 
of the proton performed in 16]. In the calculation, the quarks in the nucleon and the A have 
one unit of relative orbital angular momentum along the direction of motion (taken as the z axis) 
i.e., |AZ Z | = 1. With this we calculate the "hard amplitudes" which, in turn, are to be convoluted 
with the soft light-cone distribution amplitudes. The light-cone wave functions of the nucleon and 
delta have been classified according to their orbital angular momentum dependence in [lR 17]. 



Note that we have to add the gluon contribution to maintain color gauge invariance. We neglect 
the dynamical gluon effects proportional to the gluon field strength F a/3 which seems numerically 
suppressed in general 0]. 

We remark that a detailed calculation of the hard amplitudes for G- will involve two units of 
quark relative orbital angular momentum. Due to computational complexity it will not be pursued 
here. Nevertheless, it can be easily shown that G_ is suppressed, in the high Q 2 limit by 0(1/0;) 
relative to Go 0,0. 

Assuming our calculation for Go is relevant at Q 2 that is currently explored at Jefferson Lab, 
our result and a recent experimentally measured values of Rsm enable us to make a phenomeno- 
logical prediction of this ratio for higher values of Q 2 . We mention in this regard that the double 
logarithmic term, log 2 (Q 2 /A 2 ), will play an important role in our analysis. The importance of 
log 2 (Q 2 /A 2 ) has also been demonstrated recently in the Jefferson Lab data on the nucleon elastic 
form factor Q 2 F 2 /F 1 0. 

Our plan of representation is as follows: in section II we give the notations and definitions 
needed for the coming sections. In section III we present in some detail the calculations performed 
and comment on the obtained results. In section IV we perform a phenomenological analysis and 
make a prediction for Rsm ■ Section V will serve as a summary . 
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II. KINEMATICS AND NOTATION 



In this paper, we treat the bare delta as if it were a bound state of QCD, although in experiment 
it appears only as a resonance in certain scattering cross sections. Our analysis is relevant for 
scattering at the resonance peak where the contribution from the background vanishes by definition. 
We will not consider the so-called the dressing of the resonance due to the pion cloud which might 
be important at small Q 2 [ill. H^|. 

Consider the following scattering: 



P(P N )+ 1 *(q)^A(P A ). 



(2) 



which appears as a sub-process of the electro-production of pions. We are interested in the Go 
matrix element, in which case the photon is logitudinally-polarized and with a large virtuality 
Q 2 = —q 2 . Our analysis will be carried out in a frame in which the virtual photon 7*, the 
incoming proton and the outgoing A are collinear. The three momenta of the proton and the A 
are in opposite directions and the 7* is moving in the — z direction. The photon polarization vector 
is given by: 

■1 



e(X = 0) 



= (g d ,0,0,g u 



(3) 



with = (q°,Q,0,q 3 ). 

We will also need the light-cone wave functions of theproton and the A + (with the same charge 
as the proton). These wave functions are given by 3) 12 ] : 



|P,A = l/2) 



1 

12 

+ 



d[l][2][3]e« fc {4 1) (l,2 J 3)«J r (l 
(fc+ L ^g ) (l,2,3) + fc+ L V? ) (l,2 J 3) 
«t r (l)«} i (2)4 i (3) - ^(1)^(2)4^3)] |0>} 



4 X (2)4 T (3)- 4^(2)4^3)1 |0) 



+ 



(4) 



|A+,A 



-1/2) 



-i|d[l'][2'][3']e^{#(l',2',3') 

x «J i (l')^(2')4 T (3') +4 i (l')4 i (2')4 T (3') +d] i (l')ut i (2')4 T (3')] |0> 

+ ^IV4V,2',3') (5) 

x 4 i (io4 T (2')4 T (3') + 4^104^2')^ ( 3 + 4i(i')4t( 2 ')4 t (3')] io>} + ... 

where /cjjj = k x ±ik y , the ellipses denote other components of the wave functions which do not enter 
the following calculation, and j (J.) on the quark creation operators denotes the positive (negative) 



helicity of the quarks. The amplitudes ipp\ (rp a> ^p') have zero (one) unit of the orbital angular 
momentum projection. ^£^(1,2,3) is symmetric in 1 and 2. The integration measure is given by: 



d[l]d[2]d[3] 



dxidx2dxs d 2 ki±d 2 k2±d 2 k3± 



(2vr) £ 



x2tt5{1 -XX-X2- x 3 )(2vr) 2 ^ 2 ) (k 1± + k 2± + k 3J _) 



(6) 



and the Xi are the fractions of the proton linear momentum carried by the quarks, and ki± are 
the corresponding transverse momenta. For the A we use yi instead of Xj and k' i± instead of ki±. 
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It should be noted that since the A moves in — z direction then the orbital angular wave function 
must be k'^_ as shown in Eq. (5). The matrix element Go has two contributions: one in which 
the proton has l z = and the A carries l z = — 1 and in the second one the proton has l z = 1 
and the A carries l z = 0. For the first case, we introduce the hard amplitudes Tj(l, 2, 3, 1', 2', 3'), 
which represent three quark scattering off the external photon with two-gluon exchange. The 
index i = 1,2,3 indicates that the photon is attached to i-th quark. For the second case, the hard 
amplitudes are denoted as T/(l, 2, 3, 1', 2', 3'), in which all quark helicities are reversed. We will 
explain how to compute these hard amplitudes in the next section. 



III. LEADING-ORDER PQCD FACTORIZATION FORMULA FOR G 



In this section, we derive a leading-order pQCD factorization formula for the helicity-flip Nj* — > 
A matrix element Go- As shown in Eq. (1), this is generated by a virtual photon with longitudinal 
polarization, corresponding to the Coulomb quadrupole amplitude. 

The first step in calculating Gq is to obtain the hard scattering amplitudes of partons. For this 
we follow the guidelines of Ref. |l6|. In principle, one has to consider both three-quark scattering 
with one unit of orbital angular momentum or three-quark-one-gluon scattering without orbital 
motion. Here we consider explicitly only the former and take into account the latter through the 
requirement of color gauge invariance, leaving out the contribution with the gluon field strength 
tensor. This practice of leaving out dynamical gluon contribution corresponds to the so-called 
Wendzura-Wilczek approximation in the literature A key point in the calculation is that since 
the quark masses are negligible, the quark helicity is conserved. This is due to the vectorial nature 
of the electromagnetic coupling between the struck quark and the virtual photon and of QCD 



quark-gluon coupling 14|. When a Fock component of the light-cone wave function of the proton 
contains two quarks of the same flavor and helicity then, clearly, it contributes only when there are 
two quarks of the same flavor and helicity in the A wave function. The exchange contribution that 
results due to the presence of the same two quarks has to be taken into account properly through 
the second quantized calculation. 

Using the light-cone wave functions of the proton and the A in Eqs.(4) and (5) and keeping 
terms that are linear in the quark transverse momenta, one obtains: 

Go = -^^/d[l][2][3]^ 

x [T 2 (l, 2, 3, 1', 2', 3') + T 2 (l, 2, 3, 3', 2', l') - T 3 (l, 2, 3, l', 2', 3') - T 3 (l, 2, 3, 3', 2', l')] 

3', 2')(A;+ X ^ 3) + A;^^ 4) )(2, 1, 3) 
x [T 2 (l, 2, 3, 1', 2', 3') + 7^(3, 2, 1, l', 2', 3') - T 3 '(l, 2, 3, l', 2', 3') - 2*(3, 2, 1, l', 2', 3')] }(7) 

where Tj and T[ are the three-quark scattering amplitudes introduced in the previous section. The 
color wave function of the three quarks is normalized to unity. 

The perturbative Feynman diagrams for Tj are given in Fig. 1. There are fourteen such diagrams, 
taking into account the relative position of the photon vertex and the gluon interactions. In each 
diagram we have three collinear incoming and outgoing quarks exchanging two gluons. The photon 
could interact with any one of the quarks. With the spin-isospin structure of the proton and the A 
state functions, we assume that the first and third quarks have positive helicities while the second 
quark has negative one. 
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FIG. 1: The leading pQCD diagrams contributing to the nucleon to A transition amplitudes. The mirror 
diagrams should be added. 



To find the amplitudes Tj, we let the outgoing quarks carry orbital angular momenta k' i± and the 
incoming quarks have zero orbital angular momenta, i.e., k{± = 0. From each one of the Feynman 
diagrams we write the amplitude following the usual QED and QCD Feynman rules. Since the 
calculation is performed in the collinear frame in which the particles are highly relativistic, we 
can set the masses of the proton and the A to zero. We expand the quark spinors in the final 
state and the quark propagators to first order in the transverse momenta. Then we collect all such 
contributions from the given diagram and sum up the results of all the Feynman diagrams. This 
yields Ti(k' i± ,Xi,yi,Q 2 ). To find T/(/cjj_, Xi, yi, Q 2 ) we set k' i± = and ki± / and follow similar 
steps. 

The results for Ti(l, 2, 3, 1', 2', 3') at order of 0(k ± ) are as follows: 



11(1,2,3, l',2',3' 



-9. 



Q 4 



k{~ ± W! {xi,x 2 ,x 3 ,yi, ?/2, 2/3) + k£ L w 2 (xi, x 2 , x 3 ,y 1 ,y 2 ,y 3 ) , (8) 



where Cb = 2/3 and the functions w\ and w 2 are: 

1 1 



wi(x 1 ,x 2 ,x 3 ,yi,y 2 ,y 3 ) = 



+ 



+ 



xiVixhs x\ylx 3 y 3 y\x 3 x 2 y 2 x 3 y 3 y\x\x 2 y 2 



w 2 (x 1 ,x 2 ,x 3 ,y 1 ,y 2 ,y 3 ) = 



x 3 yiy 2 y 3 x 2 xi 
1 1 



xiyixjyi 



+ 



x\yix 3 yl 



x 3 y 2 x 2 x 3 y% 



+ 



x\y\y 3 x 2 y 2 x 3 x 2 y 2 x\y\ ' 



(9) 
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Similarly for T 2 we have: 
T 2 (l,2,3,l',2',3') = 



Q 4 



k' X j_w {y 1 ,y 2 ,y 3 ,x 1 ,x 2 ,x 3 ) + k' 3 + ± w {y 3 ,y 2l yi 1 x 3 ,x 2 ,x 1 ) ,(10) 



the function wo is given by: 

w (x 1 ,x 2 ,x 3 ,y 1 ,y 2 ,y 3 ) 



+ 



+ 



x 2 y 2 x x x 3 y\ x{y{x 3 y 3 x\Xxy\x 3 y\ 



+ 



+ 



x 2 y 2 x\y\ x\y 3 y\x\x 3 y 3 



(11) 



To obtain T 3 (l, 2, 3, , 1', 2', 3') we interchange a;i,yi,fe£J[ and 23, 2/3,^3"^ and to obtain 
T/(l, 2, 3, 1', 2', 3') we replace k'^_ with kf L and interchange Xi and y^. The above results are 
consistent with those derived in Ref. 

To further simplify the above expressions, let us define &j_-intergrated wave functions: 

4>f{ Xl ,x 2 ,x 3 ) = 2 Jld^^iktMM) 

4>f(x u x 2 ,x 3 ) = 2 J[d 2 k ± ]k 3± - (k 2± Vp 3) +k 1 j.^ ) )(fc 2 ,fc 1 ,* 3 ), 

#(x!,x 2 ,x 3 ) = 2|[d 2 k ± ]ki ± -(k 2± ^ 3) +ki^ 4) )(fc 2 ,fci,fc 3 ) 

^(ll,X 2 ,X3) = 2 j[d 2 \i L ]^{kl,k 2 M) 

(f) ( £'(xi,x 2 ,x 3 ) = 2 J [d 2 k ± ]k 3± ■\n 1 xipA > (k 2 ,h,k 3 ), 
ip%\xi,x 2 ,x 3 ) = 2 J[d 2 k±}ki± ■k 1± ip < l\k 2 ,ki,k 3 ) 
where the integration measure is defined as: 



(12) 



k* 2 k 



d 2 ki±d 2 k 2 ±d 2 k 3 ± 
(2^)6 



<5 (2) (k 1± + k 2 ± + k 3 ±). 



In terms of these wave functions and by using the obtained expressions for the hard amplitudes 
we get the following expression for Gq: 



Gn 



r 2 



[dx 



(47r ° s)2 (2M N )AQ* 3 
{^A (2/1,2/2, J/3)0p (a?i, x 2 , x 3 ) w (yi,y 2 ,y 3 , x 1 ,x 2 ,x 3 ) + w (y 1 ,y 2 , y 3 ,x 3 ,x 2 , x%) 
-w 1 (x 3 ,x 2 ,x 1 ,y 1 ,y 2 ,y 3 ) - w 2 (x 3 , x 2 , x 1 ,y 3 ,y 2 ,y 1 ) 

+ip*A 4) (2/1,2/2, 2/3)0p (^l, x 2 , x 3 ) w (y 3 ,y 2 ,y 1 ,x 3 ,x 2 , x t ) + w {y 3 ,y 2 , y 1 ,x 1 ,x 2 ,x 3 ) 
-w 1 (x 3 ,x 2 ,x 1 ,y 3 ,y 2 ,y 1 ) - w 2 (x 3 , x 2 , x 1 ,y 1 ,y 2 ,y 3 ) 

+0p\xi,x 2 ,x 3 )(p*£ 3) (yi, y 3 , y 2 ) wq(x 3: x 2 , x\,y 3 , y 2 , yi) + w (x 3 ,x 2 ,x 1 , yi,y 2 , y 3 ) 
-wi(y3,y2,yi,x 3 ,x 2 ,x 1 ) - w 2 (y 3 ,y 2 ,y 1 ,x 1 ,x 2 ,x 3 ) 

+ipp\xi,x 2 ,x 3 )(p*£ i) (yi, Vi) w {xi,x 2 ,x 3 ,yi,y 2 ,y 3 ) + w (xi, x 2 , x 3 ,y 3 ,y 2 ,yi) 
-wi(y3,y2,yi,xi,x 2 ,x 3 ) -w 2 (y 3 ,y 2 ,y 1 ,x 3 ,x 2 ,x 1 ) j, (13) 



where [dx] = dx\dx2dx 3 5{l — x\ — X2 — X3) and a s = cf s jAm. We see from Eq.(13) that Go scales like 
l/Q 4 in the high Q 2 limit, consistent with the general power counting. To find the normalization 
of Go we need to know the light-cone distribution amplitudes defined in Eq.(12). For the proton, 
a set of such functions have been given in Ref. ^| based on conformal expansion, QCD sum rules 
and Lorentz symmetry. For the A no such work has been done explicitly. However if the A is to 
be treated as a three-quark state, we believe that, at the asymptotically large Q 2 , the light-cone 
distribution amplitudes of the A will have a yj-dependence which is identical to the Xj-dependence 
of the light-cone distribution amplitudes of the prot on [20l| . For the proton we have, for example, 
the asymptotic form of 4>^ is X1X2X3 and 4>^ is X1X2 |l3(]. If we perform the Xi and yi integration in 
Eq.(13) with such functions we get, due to end-point singularities, a double logarithmic divergence 
results. This divergence indicates the quarks with small Feynman x contribute significantly to the 
hard scattering. However, for very small x for which the parton longitudinal momentum is on the 
order of Aqcd> the hard scattering picture breaks down, and the above calculation is invalid. Then 
one has to add the so-called Feynman contribution which is relatively unimportant in the limit of 
large Q 2 . 

Therefore, these divergent integrals are regulated physically by a cut-off from below of order 
A 2 IQ 2 where A is some parameter that represents the soft energy scale (order of few hundreds of 
MeV). With this cut-off, the result of the momentum fraction integration will be a Q 2 -dependent 
term of the form, log 2 (Q 2 /A 2 ). Thus we can write: 

G = c' log 2 (Q 2 /A 2 )/Q 4 (14) 

where d is a numerical factor that depends on the explicit expressions of the light-cone distribution 
amplitudes. With this form of Go and with the fact that G_ is of order 1/Q 2 relative to G+ we 
will give in the next section a phenomenological prediction of the ratio Rsm in the high Q 2 limit. 

IV. PHENOMENOLOGY 

In the electro-production of the A resonance, the multipoles Ml, E2 and G2 of the exchanged 
virtual photon are the only ones that contribute 21, Ea . Esl . E3 |. One of the experimentally 
extracted quantities related to this process is the ratio Rsm = G2/M1. In order to express this 
ratio in terms of the helicity matrix elements G\ we introduce the resonance "helicity amplitudes" 
A1/2, A3 J 2 and S1/2 |2a |. These amplitudes are computed in the rest frame of the A and the 
sub-indices refer to the helicity of the A. The scalar amplitude S1/2 is relevant only for virtual 
photons. These amplitudes are related to the helicity matrix elements through the relations [2^ : 

A 1/2 = r/G+ ; A 3/2 = i]G- ; 51/2 = r ^ G ^ ( 15 ) 



where rj is some kinematical factor and | q| is the photon three-momentum in the rest frame of the 

A given b; 
have |27|: 



Ml = -~A/2 " ^3/2 (16) 
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FIG. 2: A phenomenological prediction for the ratio Rsm- 

From the results of the previous section we see that A 3 / 2 is of order l/Q 2 relative to -Ai/2 in the 
high Q 2 limit and thus it will be neglected. With this approximation the ratio Rsm becomes: 



RsM = ~QG^ 

Substituting in the last equation the expression for |q| and using Eq.(14) we get: 

' '2.4 + Q 2 \ 2 ~~ log 2 (Q 2 /A 2 ) 



(17) 



Rsm 



2.464 



Q 2 



(18) 



The pre-factor c is determined by taking A = .25 GeV and by fitting Eq. (18) with the recently 
obtained experimental values of Rsm |lCI | : 



Rsm 
Rsm 



The result is: 



Rsm 



= -.112 ± .013 at Q 2 = 2.8 GeV 2 , 
= -.148 ± .013 at Q 2 = 4.0 GeV 2 , 

.013 V-.88 + (2.4 + Q' 2 )' 1 log 2 (16Q 2 
Q 2 



(19) 



Rsm is displayed in Figure 2. From that figure it is clear that Rsm has a slow variation at the 
presently-accessible momentum transfer Q 2 at JLab which is about 6 ~ 7 GeV 2 . We remark here 
that for the E2-to-Ml ratio Rem, two units of orbital angular momentum are need to induce the 
transition, and higher powers of the logarithmic correction may arise when one tries to calculate 
the matrix element G_ . 
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V. SUMMARY 



A perturbative QCD calculation of the helicity flip matrix element of the electromagnetic 
N —* A transition has been given. We showed that the transverse momenta of the quarks in 
the proton and A is essential to obtain a finite result. An explicit calculation of the hard ampli- 
tude Ti(xi, yi, k' i± , Q 2 ) and T((xi,yi, k{±, Q 2 ) has been given and the techniques of the calculation 
were outlined in some detail. The essential steps of the calculation are to draw the relevant per- 
turbative Feynman diagrams and then to compute the contribution from each such diagram by 
expanding the quark spinors and propagators to first power of the transverse momenta. The hard 
amplitudes are then to be convoluted with light-cone amplitudes. Our pQCD results for the scaling 
behavior of G+, Gq and G_ confirm earlier scaling prediction |7j. In addition we outlined how one 
obtains the double logarithmic correction. Based on our result for G\ we gave a phenomenological 
prediction for the ratio Rsm- Roughly speaking, I-R^mI will be of order of 20% at Q 2 of order 10 
GeV 2 . 

We thank A.V. Belitsky, C. Carlson, P. Stoler, and F. Yuan for a number of useful discussions. 
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